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The Rayleigh- Taylor (RT) instability develops and leads to turbulence when a heavy 
fluid falls under the action of gravity through a light one. We consider this phenomenon 
accompanied by a reactive transformation between the fluids, and study with Direct Nu- 
merical Simulations (DNS) how the reaction (flame) affects the turbulent mixing in the 
Boussinesq approximation. We discuss "slow" reactions where the characteristic reaction 
time exceeds the temporal scale of the RT instability, r 3> iinst ■ In the early turbulent 
stage, ii ns t < t < t, effects of the flame are distributed over a maturing mixing zone, 
whose development is weakly influenced by the reaction. At t > r, the fully mixed zone 
transforms into a conglomerate of pure-fluid patches of sizes proportional to the mix- 
ing zone width. In this "stirred flame" regime, temperature fluctuations are consumed 
by reactions in the regions separating the pure-fluid patches. This DNS-based qualita- 
tive description is followed by a phenomenology suggesting that thin turbulent flame is 
of a single-fractal character, and thus distribution of the temperature field is strongly 
intermittent. 
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1. Introduction 



The Rayleigh- Taylor (RT) instability (Rayleigh 1883 Taylor 1950 Chandrasekhar 



1961 1 occurs in many natural and man-made flows. Under constant acceleration, pro- 



duced e.g. by gravity, the instability develops into the RT turbulence (Duff et al. 1962 
Sharp 19841, characterized by a mixing zone which continuously grows. Many recent ex- 



perimental ( Dalziel e^"aT1|1999 Wilson fc Andrews||2002 ) , numerical ( Young e£"oT1|2001 



2006 



Vladimirova & Chertkov 2008) and theoretical efforts (Chertkov 2003 Chertkov 



et al. 



1941 



Dalziel et al. 1999 Cook & Dimotakis 2001 Ristorcclli & Clark 2004 Cabot & Cook 



20051 - later ones in the spirit of the turbulence phenomenology (Kolmogorov 



Obukhov 1949 Corrsin 1951 Frisch 1995 1 — have been devoted to analysis of 



turbulence inside the mixing zone. The mixing fluids can also be involved in chemical or 
nuclear reactions, e.g. flame. A flame front propagating against the acceleration is modi- 
fied by the RT instability, which leads to a reactive turbulence in the mixing zone. Studies 
of the externally stirred turbulent flames belong to the general field of turbulent com- 



bustion pamkoler||1940| |Williams|[l985l |Peters||2000"| |Poinsot fc Veynante|[2005l |Borghi 



1988 Kerstein 2001 1 and have many important applications. The buoyancy driven, re- 



active turbulence is believed to be the dominant mechanism for thermonuclear burning 



in type-la supernovae (Khokhlov 1995 Gamezo et al. 2003 Zingale et al. 20051. The 



interplay of buoyancy, reactive transformation and turbulence also plays a crucial role 



in studies of fusion (Freeman et al 



and fires (Cetegen & Kasper 1996 



19771 and large-scale combustion, such as furnaces 



Tieszen 20011. Turbulence, generated by buoyancy, 
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facilitates mixing, and thus counteracts the dynamical separation of fluids (phases) due 
to reaction. To explain this competition between separation and mixing is the main 



challenge emerging in the description of reactive flows (Kerstein 2001), which also ap- 



plies in the broader context of multi-phase fluid mechanics, e.g. the RT turbulence of an 
immiscible mixture controlled by surface tension ( jChertkov et aZ.|2005 |. 

To clarify the nature of the competition between separation and mixing, let us consider 
a typical reactive RT setting. Initially, the heavy, cold reactant is placed on top of the 
light, hot product. As in the non-reacting RT instability, the mixing zone develops. Since 
locally the reaction occurs only in the mixed fluid, turbulence inside the mixing zone 
enhances the cumulative reaction rate. The effect of reaction, however, is the opposite - 
the reaction consumes mixed material, and potentially limits the growth of the mixing 
zone. While diffusion and turbulence mix the reactant with the product, the reaction 
converts the mixed fluid into the product, thus causing an upward shift of the mixing 
zone as a whole. 

Existing studies of reactive RT setting have mainly focused on the initial stages of 
instabil ity (|Lima et aL|2006 ), or on the instability restricted by container walls in exper- 
iments (jBockmann fc Mullcr 20041 and by domain sizes in simulations (Khokhlov||1995 



Zingale et al.\2005 Vladimirova fc Rosner|2005 ) . The transverse (to direction of gravity) 
confinement results in a constant (on average) cumulative reaction speed, presumably 
dependent on the domain size (Khokhlov 1995[ ). In this system fluctuations, e.g. in ve- 
locity and density fields, do not grow with time. It is still not clear whether reaction can 
stabilize the unconfincd RT instability. One suspects that the answer might depend on 
the relation between the timescale of the RT instability and the reaction timescale. 

In this paper, we study RT turbulence unconfined by walls in the presence of reactions 
slower than RT instability. Armed with DNS and building a unifying phenomenological 
description, we analyze the interplay of turbulent mixing and reaction-mediated sepa- 
ration. Specifically, we focus on predicting macroscopic and microscopic features of the 
mixing zone, such as the upward transport of the mixing zone as a whole and statistics 
of density fluctuations within the mixing zone. 

We examine the interacting effects of buoyancy and reaction in the simplest, however 
physically relevant setting, the Boussinesq approximation, where the variations in the 
fluid density are small and are related linearly to the temperature contrast. We assume 
that the pressure increase due to reaction is negligible, and that the Mach number is 
small and, consequently, the fluid is incompressible. Thus, our system can be modeled 
by the incompressible Navier-Stokes equation coupled to the advection-reaction-diffusion 
equation governing the temperature evolution. The set of equations is written as 



d t v + (vV)v 
d t 6 + (vV)9-- 



-Vp + vV 2 v - 2Ag6, 
~ 1 R(9) + kV 2 9, 



Wv = 0, 



(1.1) 
(1.2) 



where v is the flow velocity, p is the pressure, and 9 represents temperature variations 
and/or the chemical composition of the fluid. It is convenient to choose 9 equal to zero 
in the cold (heavy) phase and to unity in the hot (light) phase. Often, 9 is referred to 
as the reaction progress variable. The reaction and production of heat is represented by 



the first term on the right hand side of Eq. (1.2); we refer to R{9) as the reaction rate 



function, and to r as the characteristic reaction time, or reaction time scale. We consider 



a single step reaction (see |Poinsot fc Veynante 2005 1 , where the reaction rate is zero in 
the pure phases, i.e. R(0) = R(l) = and R = O(l) for < 9 < 1. The buoyancy force, 



i.e. the last term on the right hand side of Eq. (1.1 1, is directed along the gravitational 
acceleration g. The factor A = (p± — p-i)j{p\ + p%) <C 1 is the Atwood number, where 
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p\ and P2 are densities of the reactant and the product, respectively. We assume that 



the kinematic viscosity v in Eq. (1.1) and the diffusion coefficient k in Eq. (1.2) are 



comparable, and thus the Prandtl number is of the order of unity. 

We compare our simulations with the benchmark non-reacting case, corresponding to 
r = +oo in Eq. (1.2 1. Phenomenological theory introduced by Chertkov et al. (2003 



2005 1 and hereafter called the "phenomenology" suggests that in the non-reacting RT 



turbulence, velocity and density fluctuations are driven by buoyancy at the scales L v 
and Lg, respectively, both comparable to the mixing zone width, H. These fluctuations 
generate cascades (of energy and temperature) towards small scales. The cascades ter- 
minate at the viscous scale, rj, and (thermal) diffusive scale, r^, respectively (Chertkov 
2003 Chertkov et al.\ 20051 . The scales n and r^ are of the same order provided the 



Prandtl number is of order unity. Previously conducted simulations of the non-reactive 
RT turbulence in the Boussinesq regime ( Vladimirova fc Chertkov|2008 1 show that while 
the mixing zone width grows in time, r\ and r^ decrease with time. It was also found 
that the correlation radii of both the velocity and temperature fluctuations, L v and Lg 
grow in time as 0.05-0.17J . One special focus was on resolving the internal structure of 
the mixing zone more systematically than in the previous studies. It was shown that 
spatial correlations do not depend much on the vertical position within the mixing zone. 
In particular, this means that the relevant values of the major scales characterizing a 
snapshot of the RT turbulence depend only weakly on the height of the horizontal slice 
within the mixing zone. 



2. Preliminary Considerations 

In this paper we discuss the regime of relatively slow reaction: the time scale of the 
non-reacting RT instability, imst, and the characteristic reaction time, r, are assumed to 
be well separated, r ^> ti ns t- 

At the early stage, t < r, development of the mixing zone and the range of turbulent 
scales inside the mixing zone are weakly influenced by the reaction. Nevertheless a cu- 
mulative effect of the reaction can be seen in the overall shift of the mixing zone, as a 
whole, from the product side to the reactant side. To estimate the shift of the mixing 
zone, Zf, let us integrate Eq. (1.2 1 along z. For the first term on the left hand side one 



arrives at dtZf, while the only other nonzero contribution is associated with the reaction 
term R/t on the right hand side. Since in this regime the heat production is determined 
mainly by the mixing on the scale of the whole mixing layer, the reactive contribution is 
estimated simply as H/t. Thus, dtZf ~ H/t. Even though the overall shift of the mixing 
zone at t < t is smaller than H , Zf grows faster than H, so that the two scales become 
comparable at t ~ r. From this time onward the reaction effects become more prominent. 

The reaction becomes more effective when the turbulent turnover time at H, estimated 
as t, exceeds the characteristic reaction time r. In this regime the reactant entrained into 
the mixing zone is burned out completely, and the overall shift of the mixing zone, Zf, 
should be of order H . This transition in the mixing zone evolution is also accompanied 
by a qualitative modification of the temperature fluctuations in the interior of the mixing 
zone. In contrast, the velocity distribution at all the scales is expected to show only weak 
sensitivity to the reaction even at t > r. 

The two stages observed for t < r and t > r will be called the mixed stage and the 



segregated stage, respectively. In the field of turbulent combustion (see Williams 1985 Pe- 
ters|2000 l, the mixed regime is called "well-stirred reactor" (|Poinsot fc Veynante|2005 1 , 
"thickened flame" ( Borghi|[l98 8) or "broken reaction zone" ( Peters||2000 ) while the seg- 
regated regime is the regime of "distributed reaction zones" (Poinsot fc Veynante||2005 1 , 
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Figure 1. Schematic representation of the system evolution in log-log coordinates: the width 
of the RT mixing layer H, the large scale fluctuations of velocity and temperature L v and Lg, 
the shift of the mixing layer Zf due to reaction, the turbulent flame thickness Lf, and the 
dissipation scales 77 and r^. The stages of the early RT instability development and transition 
to turbulence are shown shaded, r is the characteristic reaction time, separating well-mixed 
and segregated regimes. In the mixed regime, the turbulent flame thickness corresponds to the 
typical width of interfaces, separating pure phase domains of size Lg. The temporal scalings are 
phenomenological propositions discussed in Sections [2] HI 



"wrinkled-thickened flame" (Borghi 1988) or "thin reaction zone" ( Peters|2000 l. Transi- 
tion between the regimes occurring at t ss r can be interpreted as crossing the Da = 1 



line on the turbulent combustion diagram (see Poinsot & Veynante 2005). Indeed, the 
Damkoler number ( Damkoler||1940 ) , Da, is defined as the ratio of the integral time scale 
to the reaction time scale, Da = t/r, which grows linearly in time for the RT turbu- 
lence. The other popular dimensionless number in turbulent combustion, the turbulent 
Karlovitz number, Ka, describes transition between the regime of "distributed reaction 
zones" / "wrinkled-thickened flame" /"thin reaction zone" and the "flamelet" regime. The 
Karlovitz number is defined as the ratio of the chemical time scale to the Kolmogorov 
time scale. For the slowly reacting RT turbulence studied here, the "flamelet" regime 
does not apply, i.e. Ka is always much larger than unity and it increases with time. The 
laminar flame regime is also not observed in the systems considered here; consequently 
the laminar flame thickness /iiam is an irrelevant scale. 

At t > t the main effect of the reaction on the mixing zone is in creation of domains 
of pure phases (9 — and 9 — 1). In these domains the fluctuations in the reaction term 
of Eq. \l.l\ dominate fluctuations in the respective advection term, leading to complete 
suppression (burning out) of temperature fluctuations. The domains are separated by 
relatively thin (compared to domain size) interfaces where burning occurs. The interface 
width, also referred to as the turbulent flame thickness, Lf, (that is the scale where the 
turnover time of velocity fluctuations is comparable to r) becomes much smaller than L v 
at t ^> t due to the dominance of reaction. 

Formation of the conglomerate of pure reactant and product domains should result 
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A 


1600.00 


160.00 


0.10 


B 


177.78 


53.33 


0.30 


C 


64.00 


32.00 


0.50 


D 


40.32 


25.40 


0.63 


E 


16.00 


16.00 


1.00 



Table 1. Laminar flame parameter for simulations discussed in the text: reaction time r, 
laminar flame thickness /iiam = 4^ r KT, and laminar flame speed, si am = 4-^/k/t, for KPP 

reaction. 



in a qualitative transformation of the single-point statistics of the temperature field 
inside the mixing zone. A single peak distribution centered around 8—1/2 transforms 
into a distribution with two peaks, related to the emergence of pure phase domains 
corresponding to 9 = and 9—1. The qualitative differences in the projected stages are 
summarized in Fig. [I] where the temporal behavior of different characteristics is shown 
schematically. 



3. Numerical Results 

The computational technique used in this work is similar to that described in Vladimirova 
& Chertkov (2008). In our simulations, we use the popular Kolmogorov-Petrovskii- 
Piskunov (KPP) model (|Fisher||T937] |Kolmogorov et a?.|[l937|, where R = 40(1 - 6) 



with < 9 < 1. The KPP model is not the only possible choice of the function R{9). 
However, compared to other models, the KPP reaction provides a relatively thick lam- 
inar front which can be simulated at coarser resolution and makes computations more 
affordable. As before, we restrict ourself to the case of Pr = 1, and chose the units where 
v = k = 1 and 2Ag = 1 . The laminar flame parameters for five simulation sets discussed 
below are shown in Table Q] 

We solved t he equations ( 1.1|1.2 1 using the spectral element code developed by Tufo & 



Fischer (2001 1. We use elements of size 30 with 12 collocation points in each direction. 



The size of our computational domain is 960 x 960 x 1440 physical units, or 384 x 384 x 576 
collocation points. We stop our simulation at time t — 128, when the width of the 
mixing layer approaches the size of the computational domain. The boundary conditions 
are periodic in the horizontal (x,y) directions and no-slip in vertical (z) direction. The 
initial conditions, taken at t — 0, include a quiescent velocity and a slightly perturbed 
interface between the domains, determined by 9 = 9o(z + S). The function 9q(z) = 0.5[1 + 
tanh(0.4z)] describes the density profile across the interface and S(x,y) is a perturbation 
having the spectrum with modes 18 < n < 48 with spectral index 0. Here the spectral 



index refers to the exponent of the wave-number (see Ramaprabhu et al. 2005 1 , and it 
describes the shapes of the perturbation spectra. 

The least processed results of numerical simulations, snapshots of the temperature field 
9, are shown in Fig. [2] they are quite informative for our intuition. In case A (r = 1600) 
the structure of the mixing layer is practically indistinguishable from the non-reacting 
case (see Vladimirova & Chertkov 20081, at least up to t — 128. In case C (r = 64) 
the reaction starts to influence the temperature distribution somewhere between t = 64 
and t = 96, where we observe emergence of regions of pure phases. In case E (r = 16) 
reactants and products are well separated by relatively thin interfaces at t = 32. Thus, 
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t = 32 t = 64 t = 96 t = 128 





Figure 2. Vertical and horizontal (z — Zf, center of the mixing zone) slices of temperature for 

three different reaction strength. 
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a transition from the mixed stage to the segregated stage occurs at t « r, in accordance 
with the preliminary discussion of Section [21 

In Fig. [2] we observe that the width of the mixing layer, defined as the vertical distance 
between bubbles and spikes, is approximately the same in all cases. Measurements of the 
half-width of the mixing layer, defined as H = J R(9) dz (where the overbar indicates 
averaging in a horizontal plane), confirm that H(t) is weakly dependent on the reaction 
rate (see Fig. [3k). In Fig. pb we plot the root mean square (RMS) of the vertical velocity 
and see that in all five cases the RMS velocity is practically the same. Moreover, all five 
cases exhibit similar internal structure of the velocity field inside the mixing layer, as we 
see from the comparison of the temperature and velocity correlation lengths, Lg and L v , 
measured in the central slice of the mixing zone (Fig.pp). Following our earlier approach 
(Vladimirova fc Chertkov|2008 1 we define the velocity correlation length, L Vl as the half- 



width of the the normalized two-point pair correlation function of velocity, f(0)/f(L v ) = 
2, with f(r) = (v z (ri)v z (ri +r))/(v z ). The temperature correlation function is defined 
in a similar way, but based on the temperature fluctuations, 9 — 9. Both temperature 
and velocity correlation lengths differ by 20-30% compared to each other for r ranging 
from 16 to 1600. Similarly to what was reported in Vladimirova & Chertkov (2008[) for 



the non-reacting case, Lg/H sa L v /H rs 0.1. The viscous length 77 = (i/^/(5|Vw| 2 )) 1//4 
decreases slowly with time in all cases, with less than 10% variation between r = 16 and 
t = 1600. The main conclusion here is that the velocity field is essentially insensitive to 
the reaction in the slow reaction regime. 

Next we look at the shift of the mixing layer Zf — J 9 dz, introduced in Section^ We 
compute the cumulative reaction, or bulk burning rate, if — J R(9)dz, and compare it 
to the growth rate of the mixing layer, H. Both quantities are shown as functions of time 
in Fig. [3]l. We see that the bulk burning rate is greater for higher reaction rates. The 
time when H and if become comparable is approximately equal to r. 

It is not surprising that the bulk burning rate increases with decrease in r. However, 
in addition to this obvious effect, we need to take into account the reaction efficiency, 
which depends on how well the fluids are mixed within the mixing layer. In Fig. [3J5 
we show reaction efficiency of the whole mixing layer, Zfr/H, and compare it to the 
reaction efficiency in the middle slice, R(zf). Good agreement observed between the two 
characteristics suggests that the fraction of well-mixed fluid within the mixing zone does 
not vary significantly with the vertical position. In fact, we had also made this observation 



earlier in our study of the non-reacting case (Vladimirova & Chertkov 2008), where the 



mixing efficiency stayed at the constant level of « 0.8 across the whole layer. In addition, 
Fig. [3^ shows the fraction of well-mixed fluid, A mix , defined as a relative area of the 
slice z = Zf occupied by fluid with 0.25 < 9 < 0.75. We see that A m i x behaves similar to 
R(zf); some differences reflect an arbitrary choice of the selected limits. The area fraction 
A m ; x , as well as the average reaction rate in the middle slice and the reaction efficiency 
characterize the transition from the mixed stage to the separated stage. In the mixed 
stage A mix is of order unity, and in the segregated stage it is a decreasing function of 
time. Numerical results shown in Fig. ^) confirm a decrease of A m j x observed for t> r. 

Since the reaction at the central, z = Zf, slice describes the bulk burning rate so well, 
we pay particular attention to the temperature distribution within the slice, shown in 
Fig. El In case A the probability distribution function (PDF) of temperature evolves in 
time almost as if there were no reaction at all. It starts from a single peak distribution 
determined by the initial conditions (slightly perturbed diffused interface). Then, during 
the linear stage of the RT instability, t < 16, the PDF transforms into the two-peak shape 
and becomes single-peaked again at t sa 48. The transformation from one to two peaks 
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128 



128 



Figure 3. (a) half-width of mixing layer H; (b) RMS of vertical velocity at Zf, (c) pumping scale 
L v (dashed lines) and scale Lg (solid lines) measured at Zf ; (d) bulk burning rate Zf (solid lines) 
compared to growth rate of RT instability H (dashed dots); (e) burning efficiency measured as 
ratio i/r/H (solid lines), the average reaction rate R(zf) in the middle slice (dashed lines), and 
well-mixed fraction of the slice area, A m ; x (dotted lines); (f) turbulent flame thickness in the 
middle slice obtained using "expanding circle" technique. Arrows in (d)-(f) show the reaction 
time. 
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corresponds to the transition to the non-linear regime of the RT instability, associated 
with secondary Kelvin- Hclmholtz type shear instability and formation of RT mushrooms, 
while transition from two peaks to one corresponds to the destruction of RT mushrooms 
and formation of the turbulent mixed zone. When reaction is stronger, as in cases B and 
C, we observe yet another transformation in the PDF, namely the appearance of a narrow 
peak around 9 = and decrease in the middle of the 9 range at t « r. The transition 
to the two-peaked shape at t rs t is of interest to us since it indicates transition to the 
segregated regime. (In cases D and E, this 1-2-1-2 peak pattern is not noticeable because 
t is of the order of, or shorter than, the transition time from linear to non-linear RTI.) 
Quantitatively the process of the transition to the segregated regime can be described 
by integrating the PDF function over some interval in the vicinity of 9 = 1/2. Selecting 
the interval 0.25 < 9 < 0.75, we obtain A m j x , described earlier. 

The well- mixed fraction, A m ; x , is related to the turbulent flame width, Lf, or the width 
of the layer separating the pure-phase domains. Aiming to study the geometry of the 
turbulent flame we measure Lf in the following way: for each point with 1/4 < 9 < 3/4 
we find the largest radius of the circle around this point containing only points with 
1/4 < 9 < 3/4; this length, averaged and multiplied by four, represents the typical width 
of the burning region, Lf, shown in Fig.pf. We observe that (i) Lf increases with r, (ii) 
Lf ~ Lg for t < t (compare Fig. (3t and Fig. pjfc), and (iii) for t > r, Lf grows with time 
but slower than L v , so Lf/L v decreases with time. 

To summarize, we have observed in our numerics the transition at t s=s r from mixed 
to segregated regimes in variety of ways, from visual comparison of temperature fields 
and temperature PDFs to direct measurements of reaction efficiency and turbulent flame 
thickness. 



4. Phenomenology 

This Section focuses on analysis of the asymptotic (t^>r), segregated regime, suppos- 
edly characterized by extended inertial interval and scaling behaviors. The simulations 
discussed above do not run long enough to reach and accurately resolve this regime. How- 
ever, the simulations are in a qualitative agreement with the phenomenology constructed, 
thus providing a starting point for more ambitious simulations in the future. 

We first review the relevant facts from non-reacting RT turbulence, (r = oo). In this 
case the phenomenology ( Chertkov|2003 1 , supported by numerical simulations (Cabot & 



Cook 2006 Zingale et al. 2005) and some (though limited) experimental observations, 
can be summarized as follows. At t 3> ti ns t the width of the mixing zone grows as H ex t 2 . 
Velocity fluctuations within the mixing zone are described by the Kolmogorov cascade, 
whereas temperature fluctuations follow the Corrsin-Obukhov cascade. Both cascades are 
driven by the largest RT scale (the spikes and bubbles) at L v « Lg ex H . Even though this 
driving scale grows with time, turbulence at smaller scales is adiabatically adjusted to 
the growth. The adiabaticity is due to the monotonic decrease of the turbulence turnover 
time with the scale. 

As discussed above and illustrated in Fig. [I] for t < r the development in the reacting 
case proceeds as in the benchmark non-reacting case. Here 9 is well mixed within the 
mixing zone, and the single point distribution of the temperature is peaked around the 
median value of 9 = 1/2. Typical fluctuations of the velocity and temperature inside 
the mixing zone (measured in terms of the differences 5 r v and S r 9 between two points 
separated by a distance r) are described by the Kolmogorov-Corrsin-Obukhov estimates 

S r v^(L v /t)(r/L v )^ 3 , M~(r/£„) 1/3 , if L v >r>7?,r d ; (4.1) 
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Figure 4. PDF of temperature in the middle slice. Panels (a)-(e) show time evolution of 
PDFs. Panel (f) compares PDFs at t = 128 for different reaction times, r. 
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which are insensitive to the reaction. The viscous and thermal diffusion scales, r\ and r^, 
can be estimated as r\ ~ r^ ~ (ist) 3 ' 4 /L v oc i -1 ' 4 , provided the Prandtl number is of 
order unity. 

Since the velocity correlation functions are formed by the direct cascade determined by 
the non- linear term in Eq. ( 1.1 ), the driving effectively occurs at the scale L v rs Lg 



with 

the former one being not very sensitive to the reaction statistics, and thus determined 
by the same estimates as in the non-reacting case, Eq. \A.1\ . (Weak sensitivity of the 
velocity statistics to the reaction was also observed in our simulations, see Fig. [3]) In 
contrast, the 6 statistics are strongly modified at t > t when pure phase regions are 
formed separated by relatively thin interfaces (see Figs. [3] and EJ. The interface width Lf 
is estimated by balancing advection and reaction terms in Eq. (1.2 1, S r vS r 9/r ~ b r 9JT. 
One thus derives the following estimate for the turbulent flame thickness, 



L 



f 



Le{v Ll T/L 9 )W~< p */W*. 



(4.2) 



The turbulent flame thickness does grow with time, but in such a way that both Lg/Lf 
and Lf/rj also grow as time advances. (Although our numerical data do not run long 
enough to illustrate the increase of Lf with time, we did observe the increase of Lf with 
t. See Fig. |f:) 

Let us now exploit the Kolmogorov-Obukhov relation, calculating the time derivative 



of 0(ri)9(r 2 ) in accordance with Eq. (1.2 1, and then averaging the result over a slice 
perpendicular to the Z-axis for a fixed value of r = n — r 2 , where Lf <C r <C L. Making 
use of the adiabaticity and neglecting the reactive and diffusive terms in Eqs. ( 1.1|1.2 | 
one derives 

V («n) - v(r2))(0(n) - (0(ri)»(0(r 2 ) - <0(r 2 )))) - e 6 , 

where the flux term on the right-hand side originates from driving at the large-scale, 
d z (0), due to the large scale temperature gradient set by buoyancy (Shraiman & Siggia 



1994]). Using the Kolmogorov estimation (4.1 I for S r v, one arrives at 

eee-^r 2 ^ ~ {r/Hf'\ 



S 2 {r) = ((6 r 0y 



(4.3) 



The most important assumption, made while deriving Eq. (4.3), concerns neglecting the 



reactive contribution into the 9 2 balance. This assumption will be justified later in the 
Section. 

At r ^ Lf, Eq. (4.3 1 also describes the probability for two points separated by r to fall 
in two distinct domains of 8 rs and 9 s=s 1 , or vice versa. Exactly the same probabilistic 
arguments apply to a temperature structure function of any positive order, thus resulting 
in the asymptotic independence of the respective scaling exponent of the order: 

{r/H) 2 ' 3 . 



s n { r ) = (\s r e\ r 



(4.4) 



The expression implies strong intermittency of the temperature field, as the ratios S n /{S2) n ^ 2 
(H/r)( n ~ 2 ^ 3 all grow with scale at n > 2. Notice that these arguments also apply to the 
immiscible RT, of the type discussed in Chertkov et al. (2005). 

The expressions (|4.3 4.4 1 suggest that the flame interface is fractal and, moreover, 



single-fractal (as opposed to multi-fractal, see e.g. discussion of Sreenivasan et al. (1989) 
and Kerstein ( |1991| ). The fractalization of the flame is properly explained by dependence 
of the fraction, A m ; x , of the mixing zone on the turbulence flame thickness, Lf. Taking 



Eqs. 
that Zf 



(4.3 



4.4) at r — » Lf, one estimates, A m i x ~ (Lf/H) 2 ^ 3 ~ r/t. This also implies 
H, thus confirming the general conclusion made in Secti on [2[ 
Now we are ready to justify the approximation leading to Eq. (4.3 1. Let us estimate 
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contribution of the flame/reaction term, fj, = ([0(ri) — (6(ri))]R(r2)) /r, into the afore- 
mentioned 9 2 balance relation. The product is non-zero only if r 2 falls inside the interface. 
Then fj, is much less than (R)/t since 9(ri) — (9(ri)) is of order unity and it also has 
zero average. The average (R) is determined by the fraction of the well mixed region 
within the mixing zone, A m i x ~ (Lf / H) 2 / 3 . One concludes, that, indeed, /i <C 1/t ~ eg 



at r > Lf, thus justifying the assumption made above in deriving Eq. (4.3) and subse- 



quently Eq. (4.4) 



Turning to discussion of the range of scales smaller than Lf but larger than r^, one 
notices that the reaction term does not contribute to the 8 2 balance, as it is much smaller 
than the respective advection contribution. At these scales, where the direct effect of 
the reaction term is irrelevant and the Kolmogorov velocity estimates still survive, one 
expects that the Obukhov-Corrsin scaling, S r 9 ex r 1 ' 3 , applies. Matching the self-similar 
behavior with Eq. (4.4) at Lf 7 one derives the following scaling relation for the structure 



functions at Lf ^> r ^> Td 

S n (r) ~ {L f /H) 2 / 3 {r/L f r' 3 . (4.5) 

This formula allows a simple interpretation: the first term on the right hand side stands 
for the fraction of the well mixed region within the mixing zone, A m j x , where the reaction 
takes place, while the second term accounts for the Corrsin-Obukhov decay of correlations 



with the decrease in scale. The main contribution in the RHS of Eq. (4.5 1 comes from 
the interfacial domains, while contributions associated with pure phase domains, where 
9 sa or 9 « 1, is much smaller. 



5. Conclusions 

In this manuscript we analyze the reactive RT turbulence in the case of a slow flame 
realized when the typical reaction time r is larger then the time of RT instability de- 
velopment, tmst- The RT instability leads to formation of the mixing zone at t ~ ti nst . 
Further development of the mixing zone is roughly split into the early turbulent stage, 
t >t> ti ns t, characterized by well-mixed hot and cold phases, and the later stage, t > r, 
where pure state phases are segregated. In the early stage, the system is only mildly sen- 
sitive to the reaction (flame), resulting in a slight shift of the mixing zone as a whole, but 
no visible feedback of the flame on velocity distribution. The velocity fluctuations are 
still insensitive to the reaction, at t > t, however the temperature fluctuations become 
modified significantly. In this regime, burning takes place within thin interfaces separat- 
ing large patches of pure phases. Our numerical simulations confirm the basic prediction 
of the transition at t ss t. All the qualitative features expected from the transition are 
observed in the simulations. Given that our numerical resources were not sufficient for 
detailed analysis of the asymptotic, t>7", stage, we rely here on the phenomenology. 
Separation of the advective range for density fluctuations in two ranges, and formation 
of a single-fractal flame interface are our two main phenomenological predictions. This 
highly intermittent behavior of density fluctuations is in a great contrast with what 
was observed without reaction. These and other predictions of the phenomenology are 
subjects for juxtaposition with future numerical and laboratory experiments. 

Notice that the general picture of mixing at t > r discussed above is reminiscent of 



the immiscible turbulence discussed in Chertkov et al. (2005 1. In both cases large regions 



of pure phases are separated by thin interfaces. Viewed at the scales larger than the 
interface width, the interfaces are advected passively. Moreover, strong intermittency in 
the density fluctuations, reported above for the reactive case, also takes place in the 
setting of externally stirred (e.g. by gravity) immiscible turbulence. 
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Finally, our phenomenology also applies to other regimes of turbulent flames, e.g. those 
realized in combustion engines, where the turbulent flame width is positioned in between 
the integral scales of turbulence and the viscous/diffusive scales, L v , Lg 3> Lf 3> n, r^. 
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